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Abstract 

A model for lattice fermion is proposed which is, (i) free from doublers, (ii) 
hermitian, and (iii) chirally invariant. The price paid is the loss of hypercubic and 
reflection symmetries in the lattice action. Thanks to the e-prescription, correlation 
functions are free from the ill effects due to the loss of these symmetries. In weak 
coupling approximation, the U(l) vector current of a gauge theory of lattice fermion 
in this model is conserved in the continuum limit. As for the U(l) axial vector 
current, one obtains the ABJ anomaly if the continuum limit is implemented before 
the chiral limit m = 0. The anomaly disappears, as in the Wilson model, if the 
order of the two limits is reversed. 
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1 Introduction 



Lattice formulation provides a framework for quantum field theory where one can, in 
principle, address questions which are non-perturbative. Interacting gauge fields are easily 
incorporated through the link variables and the lattice spacing 'a' provides built-in gauge- 
invariant regularisation. In practical applications, however, a serious impediment is the 
problem of fermion doubling 1 . The 'naive' lattice action leads 1 to fifteen 'doublers' all 
with the same mass as the physical fermion and survive in the spectrum in the continuum 
limit. 

A 'no-go' theorem 2 states that the unwanted doublers can be evaded, but only at the 
price of some basic, desirable properties 3 of lattice fermion action. A variety of models 
have appeared in literature which get rid of the doublers, either partially or completely, 
by abandoning one or more of these basic properties. The most popular among these, 
the Wilson model, abandons chiral symmetry 1 ' 4 . The Wilson term in the lattice action of 
this model vanishes formally in the continuum limit, i.e. is an irrelevant term, and breaks 
chiral symmetry explicitly. It does the job of removing all the doublers by giving each of 
them a mass of (1/a), i.e., of the order of the cut-off. For a vector-like gauge theory, 
e.g., QCD, the Wilson model is adequate. In weak coupling perturbation theory (wept) 
it reproduces in the continuum limit the correct Ward identities, and, most importantly, 
the ABJ anomaly in the U(l) axial current 1 ' 5 . The explicit breaking of chiral symmetry, 
however, renders the Wilson model unsuitable for chiral gauge theories like the Standard 
Model. The difficulty with chiral gauge theory is not specific for the Wilson model, but 
an outstanding problem of lattice formulation and quantum field theory in general. 

It is clear that in order to be compatible with chiral gauge theories the lattice action 
for fermion should have exact chiral symmetry. Even the irrelevant terms, like the Wilson 
term, needed to remove the doublers from the spectrum, must leave chiral symmetry 
unscathed. In view of the no-go theorem, it is also clear that the chirally symmetric 
irrelevant term must pay some price. What, however, is not clear is whether the model 
will still be physically acceptable. A decisive criterion for this would be the Ward identities 
in wept. We report here the results of a search for a chirally symmetric lattice fermion 
action which satisfies this criterion. 

In sect. 2 we introduce various options for finite differences on lattice, all of which are 
candidates for derivative in continuum. The specific choice of finite difference for our 
lattice action (sect. 3) is motivated by the three principles, (i) absence of doublers, (ii) 
strict chiral symmetry, and (iii) hermiticity. The price paid 3 is the loss of hypercubic and 
reflection symmetries 6,7 . We, however, give a prescription to remedy the lack of these 
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symmetries in correlation functions. In sect. 4 we derive the Ward identities of the U(l) 
curents in wept and show that the vector current is conserved while the axial current 
has the ABJ anomaly in the continuum limit. Important steps in the derivation of these 
identities are given in the Appendix. We conclude with some remarks in sect. 5. 



2 Finite Differences and Fermion Doublers 



There are several options for finite differences on lattice, all of which coincide with deriva- 
tive in the continuum limit : 

(a) forward difference 

K = = -iP„ (2.1) 

(b) backward difference 

1 _ P W^ 

% = = ~ iP t ( 2 - 2 ) 

(c) symmetric difference 

K = 7T^— = (2-3) 

" 2a 2 y ' 

where p^s are generators of translation. The forward and the backward differences may 
be expressed in terms of the symmetric difference 

si = ^ + a«, s; = s; - a* (2.4) 

where the antisymmetric combination A£. 

A; S ±_ (2.5) 

vanishes in the continuum limit. 

The 'naive' lattice action is the result of sustituting the symmetric difference 8^ for 
the derivative <9 M in the continuum euclidean action 



Ta 



(2.6) 



X,/i 
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where \i = ae M with a unit vector in the /i-direction. Within the Brillouin zone 
—n /a < < Ti /a the inverse propagator 

(X 

derived from the naive action vanishes at k^ = and k^ = ix/a. The entire Brillouin zone, 
therefore, splits up into 16 distinct sectors corresponding to the domains 

< < f , f < \K\ < - (2.8) 

2a 2a a 

of each component of the four momenta, in each of which the fermion propagator assumes 
the standard continuum form 

G™ wl {k)- 1 = irfk ll 

with 7^(A = 1,2, ....16) unitarily equivalent to 7 M . The physical fermion resides in the 
first hypercube < \k^\ < ir/2a, all //, and corresponds to 7^ = 7^, while the 15 doublers 
live in the remaining sectors 1 ' 4 . 

Appearance of the doublers is not peculiar to the naive action, but, according to the 
no-go theorem 2 , is generic. A simple version of the theorem states 3 that the 15 doublers 
are unavoidable if the lattice inverse propagator has (i) continuity, (ii) hypercubic symme- 
try, (iii) reflection symmetry i.e., 7467^ (A;)^ = G F l (hi — k^j, and (iv) chiral symmetry 
75G^ 1 (/c)75 = —G F l (k). Models have been proposed which eliminate the doublers, either 
completely or partially, at the expense of one or more of the above properties. 

In Wilson model chiral symmetry is abandoned to achieve complete elimination of the 
doublers in the continuum spectrum 

(p-ip^a _ p ip^a\ / p-ip^a _i_ pip^a _ o\ 

—*r-) +r { * K 

= - X) 7T [$i x + A*) ( r + 7 M M^) - ij>(x) (r - 7 M )^(x + n) ( 2 - 9 ) 



The Wilson term makes a mass-like contribution of (r/a) in the inverse propagator 



c^ikr 1 = e 



sin(/c„a) r 
%ln V - (1 - cos(A^a)) 

(X (X 



(2.10) 



in the neighbourhood of k^ m 71 /a where the doublers live. This lifts the degeneracy 
between the physical fermion and the doublers, and removes the latter from the spectrum 
in the limit a = 1 ' 4 . 
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3 Chirally Invariant Model for Lattice Fermion 



The species doublers in the naive lattice action (2.6) owe their origin to the substitution 
of the symmetric difference for continuum derivative. A recipe for cure, which follows 
naturally 7 , is to use instead forward (2.1) or backward (2.2) difference. A serious fallout 
of this recipe is the loss of hermiticity of the action. The integrand in the partition 
function ceases to have the interpretation of a probability measure, and Feynman rules 
would yield terms which do not have correct reality properties in the continuum limit 4 . 

Intimately related to the hermiticity of action is the chiral structure and anomaly of 

a gauge theory. Consider the continuum Dirac operator of a non-abelian chiral gauge 
theory 

iP±(A) = - (f - i9^A) (3-1) 
with = A^. The effective actions V + (A) and T-(A) defined by 

T±(A) = - In det (ip±(A)) (3.2) 
are related to each other through hermitian conjugation 

r + (A) = r_(A)t (3.3) 

The relation (3.3) is of essence for the chiral properties of the theory 8 and implies that 
chiral anomaly is related to imaginary part of the effective actions T±(A). Eq.(3.3) is a 
direct consequence of the hermiticity of the Dirac operator p(A). 

p(A) = (i$+gA) 

( D(A) \ (3.4) 
" [Dl(A) ) 

The Weyl components D(A), D\A) 

D(A) = a II (td II + gA fl ) (3.5) 
with a M = (i,<Jk), are hermitian conjugates of each other and eq.(3.3) follows from 

F+(A) = - In det [D(A)D*(0) 

r_ (A) = -In det [,0(0)^(^)1 (3.6). 
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The naive lattice action (2.6) with gauge field interactions through link variables U l 



gnaive = _ £ < ^ [ ) ^ > ( 3 .7) 



I- 1 



obeys the hermiticity condition even for finite lattice spacing. The fact that it fails to re- 
produce the AB J anomaly is well understood 1 ' 4 in terms of cancellation of the contribution 
from the physical fermion by those from the doublers. 

The key elements for a chirally invariant formulation of lattice Dirac fermion, therefore, 
are (i) absence of doublers (ii) chiral symmetry, and (iii) hermiticity. Note that there is 
no requirement of hermiticity for the Weyl components of a Dirac operator. All these 
elements are realised if in the continuum Dirac operator we substitute forward (backward) 
and backward (forward) finite differences for the derivatives in the right-handed and the 
left-handed Weyl components respectively 

4H=) o° J (3 ' 8) 

In Dirac basis the wave operator (3.8) assumes the form 

(p-ip^a _ p ip^a\ / „-ip M a , „ip M a o\ 
Ta ) " V» ( ~a ) = ^ " ^ < 3 ' 9 ' 

The 'irrelevant term' 7 M 7 5 A^ coincides with the Wilson term in (2.9) if in the latter the 
parameter r is replaced by 7^75. The representation (3.9) has the correct hermiticity 
properties, is chirally invariant, and, as is easily verified, free from doublers. 

The price paid 3 for eliminating the doublers is the loss of hypercubic and reflection 
symmetries. Hypercubic symmetry is the remnant of SO (A) symmetry on lattice and 
is the invariance under rotation through n/2. This is lost because the antisymmetric 
combination A^ does not transform as a vector. As for reflection symmetry, 7475A4 is 
compatible but the remaining pieces in the irrelevant term are not. Note that the sum 
of squares ^(7 M 7 5 A^) 2 , does not suffer from these incompatibilities. This suggests as 

remedy to use instead the wave operator 

- 7^ M 75 = 7^' - 7^75 A« (3.10) 

where e M = ±1, with the prescription that all correlation functions are to be obtained after 
averaging 7 over e M . What e M achives is the decoupling of the finite difference, forward or 
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backward, used for the /z-th component of the derivative from those used for the other 
components. The (anti) correlation of finite differences used for the right-handed and the 
left-handed Weyl components is, however, maintained by the wave operator (3.10), so 
that the chiral structure (3.4) is preserved. The wave operator (3.10) leads to the 'free' 
fermion propagator 



sin [k^a] 



+ 7^75 



1 — cos (k^a) 



(3.11) 



Interactions with gauge fields are introduced, as usual, through the link variables U^, 
and the lattice action corresponding to the wave operator (3.10) is given by 



S F = 



J2^<^\h(^ a ut- U ^ a ) 



m 



+7*75 (2 



(3.12) 



-i)(x)Y^U^(x)i)(x + n)} - mijj(x)ip(x) 
We have introduced a mass £ m' for the fermion which vanishes in the chiral limit. 

Ward identities for the U(l) vector and axial vector currents are derived by requiring 
invariance of the partition function under the local transformations 



ip(x) 



Ja(x 



ty(s), 



i){x) -> ^(x)e- ia{x) 



and 



if>(x) -> e t/3(xh ^{x), ij)(x) -> ^{x)e i/3(x) ^ (3.13) 

respectively. The Jacobians of fermion measure in the lattice fermion action are trivial 9 for 
both the transformations (3.13). Anomalies, if any, can arise only from lattice artifacts, 
the 'irrelevant terms'. The Ward identities on lattice, therefore, are 



a 

1 

a 



where 



< J%( x ) ~ J U X -V)>=-< J l ( x ) - J l ( x ~ A*) >> 
J t( x ) = \ Up{xh^{x)^j{x + fi)+ $(x + /j,)^Ul(x)i)(x) 



(3.14) 
(3.15) 



J U ( x ) = \ [^(x)%U^x)^(x + ip(x + n)^Ul(x)il>(x) 

J^( x ) = l - [^l^U^ix + n)- ${x + tiifaUfoMx)] . (3.16) 

The right hand sides of eqs.(3.14) and (3.15), if non-zero, are to be identified as anomalies 
in U(l) vector and axial vector currents respectively 



< d^J^x) > = lim < Y > 



< d^J^(x) > = lim < X > 



= lim - < Jf. (x) — Jf. (x — a) > 



(3.17) 



(3.18) 



4 U(l) Ward Identities in WCPT 



To evaluate < X > and < Y > it is convenient to start from their operator representations 

< y >= Tr lb < x\ (G F f? + fi'Gp) \x > 
<X>=Tr<x\ (G F f? - fi e G F ) \x > 

where ffi = 7a£a-Ra with 



2 - U x e ip * a - e~ ip * a U\ 



2a 



The fermion propagator is given by 

Gp 1 =i(p- ifi e "f 5 - im) 

where 



D, = 



U x e ipxa - e~ ip ^ a U\ 
2ia 



(4.1) 
(4.2) 



(4.3) 



(4.4) 



In (4.1) 'Tr' stands for trace over 7-matrices, and, in non-abelian gauge theory, also over 
the symmetry matrices. 
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In wept the link variable U x ~ 1 + iagA x , so that for small lattice spacing 

D x = (s\ + gA x cos(p x a)) + 0(a) , . 

= c A + 0(a) ^ 

where 



sin(p A a) l-cos(p A a) 

sa = , c A = (4.6) 

a a 



The commutators in u>cp£ are given by 



[R\, R p \ = [D x , R p ] = ^ ^, 

[-D A , -Dp] = F Ap cos(p A a) cos(p p a) = / Ap 

The quadratic form of the propagator Q, defined by 

Q- 1 = (GpG^y 1 

= (p + iR e >y 5 - im) (p + iR e ^ + im) 

= ipl + ^a) + ^x P fx P - 2 l5 a Xp e p R p D x + 



(4.8) 



,2 



commutes with 75 and its trace with odd powers of 7-matrices vanish. These properties 
lead to a simpler structures for < Y > and < X > 

<Y> = Tr l5 <x\(pGR e + R e Gp)\x> 

<X> = Tr <x\(pgR e - R e gp)\x> -2iT n5 < x\R e GR e \x > 1 } 

In wept where both the lattice spacing 'a' and the gauge coupling are regarded as 
small, one can define a 'potential' V 

g- 1 = g- 1 + v 

and develop a perturbative series 5 

g = g - g vg + g vg vg + (4. 10) 

where 

<? - 1 =E(*a + <$)+™ 2 (4- 11 ) 

A 

and 

V = (ZgAxSx cos ( Px a) + g 2 A\) + ^ ]T a Xp f Xp - 2 7s ]T a Xp t p R p D x (4.12). 

A A,p A,p 
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Two factors play key role in the calculation of < X > and < Y > : (i) terms odd in R x 
drop out because of e-averaging, and (ii) in the 'physical' sector of the loop momentum 
< \k\\ < n/2a,R x is of 0(a) whereas in the 'doubler' sector 7r/2a < \k x \ < n/a it is of 
(1/a) and behaves like the mass of a Pauli-Villars regulator field 10 . One obtains (see 
Appendix) 

lim < Y >= (4.13) 
so that the U(l) vector current is conserved 

(4.14) 
In the expression for < X >, however, only the first term in (4.9) vanishes so that 



< d p J^(x) >- 



< X > = -2iTr 75 < x\R e GR e \x > 

^cfjl a (cos(k a a)) 

= 2itr(F Xp F Xp )f^-$^f 

J2(sl + c 2 x )+m- 

4 

= (i/2n±)tr(F Xp F Xp )J2(-iy 3 Cv-iIv 

u=l 

where F Xp is the dual field tensor, 'tr' stands for trace of symmetry matrices, and 

^ l+cos(fc M a) 2 



(4.15) 



<?k- 



E s a + E — - +m 



(4.16) 



_A=1 A=l 

In the continuum limit l v = n 2 /2u. Thus the model reproduces the ABJ anomaly 



lim < X > 



i/16n 2 )tr (F Xp F Xp ) 



(4.17) 



Note that in (4.9) the subscript v of I v has the meaning of the number of components 
of the loop momentum with support in the doubler sector n/2a < \k^\ < ir/a. For finite 
'm' the contribution from v = vanishes in the continuum limit. If, however, m is zero 
the latter would exactly cancel the right hand side of (4.15). For a nonvanishing anomaly 
it is, therefore, essential that the continuum limit is taken first and chiral limit m — 0, 
if necessary, afterwards. This is true also for the calculation of anomaly in the Wilson 
model 5 . 
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5 Concluding Remarks 



Our search for a model of lattice fermion which is, (i) free from doublers, (ii) hermitian, 
and (iii) chirally invariant, has paid rich dividends. Not only is the model free from 
doublers, it reproduces in wept in the continuum limit the correct Ward identities, and, 
most importantly, the ABJ anomaly in the U(l) axial vector current. In this respect the 
model measures up to the Wilson model. In a very crucial way, however, the two models 
differ. Whereas chiral symmetry is broken explicitly by the Wilson term, the irrelevant 
term in the proposed model is chirally invariant. This makes the latter a viable candidate 
for chiral lattice gauge theories. 

An important lesson of the present exercise is that the ABJ anomaly can be realised 
even with a chirally invariant term. What is important is that the physical fermion has, 
to start with, a small but non-zero mass and the chiral limit m = is taken only after the 
continuum limit. If the order of the two limits is reversed the ABJ anomaly disappears. 
The same observation, interestingly enough, is true also in the chiral symmetry breaking 
Wilson model. In the absence of the mass term, a chirally invariant action for lattice gauge 
theory breaks up into two completely decoupled Weyl components. In Weyl basis, vector 
and axial vector currents lose their individual identity and both the currents are conserved. 
By coupling the two Weyl components, the mass term brings about an asymmetry between 
the vector and the axial vector currents. This asymmetry is apparently the genesis of the 
ABJ anomaly. The ABJ anomaly, once it materialises, persists even when the parameter 
m, which triggered it, vanishes. There is a parallelism here with the phenomenon of 
ferromagnetism. A magnetic field, however weak, is needed to trigger the magnetised 
state of ferromagnet, and, once realised, the state persists even in the symmetruy limit 
corresponding to zero field. For ferromagnetism, we have to take the thermodynamic 
limit first and the symmetry limit afterwards just as for ABJ anomaly the continuum 
limit a = must precede the chiral limit. 
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Appendix 



We show in the following how the first term of (4.9), let us call it < X D >, and < Y > 
vanish in the continuum limit. The crucial ingredients are Dirac trace, e^-averaging and 
locality of the matrix elements. 

To illustrate the locality property, concentrate on a typical term 



In the physical region < < n/2a,c a ~ 0(a) and the integral vanishes in the 

continuum limit, whereas in the nonphysical or doubler sector 7r/2a < \k^\ < ir/a, the 
denominator in the square bracket is of (1/a 2 ) but the numerator is at most of order 
(1/a) unless x = y. The same argument can be continued and extended for the remaining 
factors and hence proves the locality prperty. Note the parallelism here of the role of c\ in 
doubler sector with that of the mass of the Pauli-Villars regulator field in the calculation 
of chiral anomalies 10 . 

We now illustrate (4.13). Using locality and Dirac trace 




< x\R a Q fx p Q RxD v Q Dp\x > 
d A k 

: > TTT^i < k \f\ P \y > dA y < y\QoR\DvQoVfi\ x > 




A 



< y > = trj 5 < x\pgfi e + R e gp\x > 

= -2tr 75 < x\R e a D a g vg vg \x > 




= 



Next we turn our attention to < X >. 



< X >=< X D > - < X\ >, 



where 



< X >= Tr < x\pgR e \x > 
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now, 

< X Q >=< X ol > + < X o2 >, 

with, 

< x o1 >= tr < x\^g v b g l -a Xp h P g p\x > 

and 

< x o2 >= tr < x\fi e g l -a 2p f 2p g v 5 g p\x >, 

after plugging in the perturbative expansion for Q in terms of Q and V. 

After e M -averaging, putting in R 2 /A in place of R 2 a ^R 2 = ^ R 2 \j > an d using locality, 
we finally obtain, 

< X ol >= - Tr ilbVpvPxp) < x\R 2 DpQ D v Q j Xp Q \x > 

/3,\,v,p 

Similarly, we find, < X o2 >= — < X ol >, so that < X Q >= 0, implying < X >= in 
the continuum limit. 
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